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Abstract 



o 

5^ ■ We find various asymptotic limits of universes containing two interacting fluids which interact 



and exchange energy. We study the finite-time singularities that may develop in these cosmolo- 



^v^j ' gies and obtain the asymptotic behaviours of the solutions near such blow up regimes. 

>' 

In Ref . PP , we studied the dynamical problem defined in the geometric setup of a flat FRW universe 
with scale factor a(t) containing two fluids with equations of state 

g : pi = (r-i)pi, P2 = (7-i)P2. (i) 

m 

In particular, we are interested in the behaviour of this system near finite-time singularities, in 
an effort to unravel the dynamical nature of blow up singularities and other asymptotic solutions 
that such a model can sustain asymptotically. The total energy-momentum tensor is given by 

*_ , 

j^totai) _ rp(i) _|_ j,{2) ^ an< ^ we consider an interaction of the form 

H 

s = — 

a z 

where a, f3, m, n, A, p are constants, H is the expansion rate and v° is the timelike direction of the 
fluid velocity. The evolution of this system is governed by the equations 

3H 2 = Pl +p 2 
pi + 3HTpi = -j3H m p\ + aH n p% (3) 
p 2 + 3H 1P2 = f3H m p^-aH n p^. 



- v ° = -f3H m p x 1+ aH n pZ, (2) 
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We restrict here to the case where all the exponents m, n,\,/j,= 1, the general case being currently 
under investigation (Ref. [2]). Renaming x = H, the Einstein equations ([3]) become equivalent to 
the following cubic dynamical system, 

x = y (4) 
y = -Axy - Bx 3 , (5) 

where A = a + f3 + 3 , j + 3T, B = 3(aT + /?7 + 2>T^)/2. The possible solutions of this system 'at 
infinity' correspond to the asymptotic orbits of the vector field 

f(x,y) = (y,-Axy-Bx 3 ). (6) 

The method of asymptotic splittings developed in Ref. [3] scrutinizes all possible modes that the 
vector field (jSJ) attains on approach to the finite-time singularity located at t = 0. These modes 
correspond to the inequivalent ways that (jBJ) splits as t — > 0, and are given by the following three 
distinct decompositions, 



;(!) 


= (y, 


—Axy — Bx 3 ), (all-terms-dominant case) 


(7) 


/( 2) 


= (y, 


-Axy), 


(8) 


f V) 


= (y, 


-Bx 3 ). 


(9) 



When applying the method of asymptotic splittings in our system, we find that we can have in 
total seven different families of asymptotic solutions, cf. Ref. pQ, depending on the Barrow-Clifton 
parameter of Ref. [4] given by 5 = B/A 2 . In Ref |T], we looked for Fuchsian type solutions, that 
is when the leading order of the formal asymptotic series is a power-law with rational exponents 
(in particular there is no constant leading term). In the limit 8 — > 0, the second decomposition 
/( 2 ) = (y, —Axy) leads to the unique balance, 

B? ) = [S,p] = [(2/A,-2/A),(-l,-2)], A/0, (10) 

and this in turn is used to build the self-consistent, general (two arbitrary constants, the position 
of the singularity and the shown coefficient of the asymptotic series) asymptotic Fuchsian solution 
valid around the singularity, 

x = l^ 1 + c 21 t - —c 2 2l t 3 ■ ■ ■ , (11) 

while the y expansion is obtained from the above by differentiation. 

In the limit 5 — > — oo, the third decomposition of the system has two balances that lead to 
analogous results, namely, 

B® = [(±y/2/^,^y/y=B),(-l,-2)], B<0. (12) 
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In this limit of the parameter 8, we meet a case of 'phantom fluids'. More precisely, this asymptotic 
solution exists provided the fluid parameters T,7 satisfy the condition: T^f < 0. The general series 
expansion is then given by the form 

x = ±y -^t- 1 + c 4 it 3 =F — c li y + • • • • ( 13 ) 

For the range of parameter values < 5 < 1/8 (this is the case of 'standard decay' of Ref. PQ), 

the fluids decay in proportion to their energy densities. The 'all-terms-dominant' case sustains the 

balances _ 

r / A ± VA 2 + 8B -At AVA 2 + 8B\ 

,(-1,-2) 



°1,2 



25 ' 25 
and these lead to the the asymptotic forms 



(14) 
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x = —t' 1 + c n + —(? n t + • • • , for and 8 = 1/9, (15) 
x = ^-t -1 , for 8 = 1/8. (16) 

The most dominant term of the asymptotic solutions (|lip - (|16p on approach to the singularity 
describes either a collapse (big bang type), or a big rip singularity. 

There are also other families of solutions. The first decomposition of the system reveals a very 
complicated singularity for the case 8 > 1/8. In this case, the system no longer admits a solution 
in a form of a Fuchsian series, but instead there is a new family of solutions given by the form, 

* = i±^ t -. (17) 

The scale factor here is an oscillatory function of time, increasing or decreasing exponentially, 
depending on the value of A. There is another interesting case describing anti-decaying fluids when 
8 — > oo, which in terms of energy transfer means that p2 decays while p\ gains energy instead. From 
the third asymptotic decomposition of the system and when B > 0, this oscillatory behaviour of 
solution is given by 

x = TiVVB t- 1 + c 4 it 3 T * ^2 VVBc 2 Al t 7 + ■■■ . (18) 

Lastly, regarding the all-terms-dominant decomposition, there is a dynamical situation stemming 
from the first dominant balance asymptotically, that is for 8 = 1/9 there is a valid series expansion 
given by, 

6 

x = -r 1 + cut- 2 + c 2 it" 3 + • • • , (19) 
and describes a universe driven away from any finite-time singularities. 



3 



It is interesting that many of these solutions exist in a generalized framework when there is 
curvature and more general interactions, and this more general problem is currently under inves- 
tigation in Ref. [2j. There appear to be new regimes that exist only under the joined asymptotic 
influence of curvature and interaction between the fluids. 
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